Scaling of the Lyapunov exponent in type-Ill intermittent chaos 
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The scaling behaviour of the Lyapunov exponent near the transition to chaos via type-Ill in- 
tcrmittency is determined for a generic map. A critical exponent expressing the scaling of the 
Lyapunov exponent as a function of both, the reinjection probability and the nonlinearity of the 
map is calculated. It is found that the critical exponent varies on the interval < /3 < 1. This 
contrasts with earlier predictions for the scaling behaviour of the Lyapunov exponent in type-Ill 
intermittency. 
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1. Introduction 



Intermittency is a common scenario for the transition 
to chaos in nonlinear dynamical systems. Intermittent 
chaos is characterised by the display of long sequences of 
periodiclike behaviour, called the laminar phases, inter- 
rupted by comparatively short chaotic bursts. The phe- 
nomenon has been extensively studied since the original 
work of Pomeau and Manneville [l[ classifying type-I, - 
II, and -III instabilities when the Floquet multipliers of 
the local Poincare map associated to the system crosses 
the unit circle. Type-I intermittency occurs by a tan- 
gent bifurcation when the Floquet 's multiplier for the 
Poincare map crosses the circle of unitary norm in the 
complex plane through +1; typc-II intermittency is due 
to a Hopf's bifurcation which appears as two complex 
eigenvalues of the Floquet's matrix cross the unitary cir- 
cle off the real axis; and type-Ill intermittency is associ- 
ated to an inverse period doubling bifurcation whose Flo- 
quet's multiplier is —1 2]. Although other mechanisms 
may occur leading to intermittency, these three cases are 
the most simple and the most frequently encountered in 
low-dimensional systems. 

Many experimental evidences for these types of inter- 
mittency have appeared in the literature. In particular, 
type-Ill intermittency has been found in lasers elec- 
tronic nonlinear devices @, H, @], biological tissues 0, 
and epilepsy The statistical signature of intermit- 
tency is usually given by the scaling relations describing 
the dependence of the average length (denoted by (I)) of 
the laminar phases with a control parameter (denoted by 
e) that measures the distance from the bifurcation point. 
Chaotic dynamics is characterised by the positive sign of 
the largest Lyapunov exponent (denoted by A), although 
this quantity is in general more difficult to measure from 



experimental data than statistical variables such as the 
average laminar length. Pomeau and Manneville [l[ as- 
sumed a uniform reinjection probability into the laminar 
phase and calculated the scaling behaviour of both the 
average laminar length and the Lyapunov exponent for 
those three types of intermittent chaos. Specifically, for 
type-III intermittency Pomeau and Manneville predicted 
the relations {1} oc e _1 , and A oc e 1 / 2 , both when e — > 0. 



However, theoretical and numerical studies 
as well as some recent experiments 
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16j . have shown deviations from the Pomeau- 
Manneville's prediction for the scaling behaviour of the 
average laminar length for type-III intermittency; they 
found (I) oc eT v , with 1/2 < v < 1. These discrepancies 
are attributed to the presence of mechanisms leading to 
nonuniform probability reinjections. On the other hand, 
the scaling behaviour of the Lyapunov exponent in type- 
III intermittency has rarely been explored. 

The aim of this paper is to clarify the effects of both 
the reinjection probability and the nonlinearity on the 
scaling properties of the Lyapunov exponent at the tran- 
sition to chaos via type-III intermittency by using a sim- 
ple one-dimensional map. We show that the behaviour of 
the Lyapunov exponent near this transition exhibits ap- 
preciable deviations from the prediction of Pomeau and 



Manneville. We find that A 



e' 3 , for 



0, where 



the critical exponent (3 varies continuously in the inter- 
val < (3 < 1 as a function of the map parameters. 



2. Scaling Behaviour of the Lyapunov Exponent 

The return map is a fundamental information that can 
be obtained experimentally in a chaotic system and in- 
evitable contains some reinjection process for type-III in- 
termittency. Therefore, a reinjection mechanism should 
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FIG. 1: Generic map Eq. ((1} for producing type-Ill intermit- 
tency, with z = 2, m — 2 and e = 0.002. 



be taken into account in the analysis of characteristic 
quantities for intermittency such as the mean laminar 
length or the Lyapunov exponent. We consider the fol- 
lowing one-dimensional map [Ic| 



x n+ i = f(x n ) = 



(1 + e)x n + x z n , 
{l-x n ) m , if x 



if x n < x* 

> x* 



(1) 



where x n € [0,1], z > 2, m > 1 and a;* is the solution of 

(1 + e)x* + x* z = 1 . (2) 

The parameter e measures the distance from the bi- 
furcation point. For e > 0, the map Eq. |T]) generates 
type-Ill intermittency near the origin which becomes an 
unstable fixed point. The parameter m controls the rein- 
jection probability into the laminar region of the map, 
with m — 1 corresponding to uniform reinjection. The 
reinjection probability becomes more densely localised 
near the origin as m increases. Figure [1] shows the map 
Eq. p. 

The Lyapunov exponent is calculated as 



lim 



1 

N 



N-l 

E 

n=0 



log|/'(x„)|. 



(3) 



We have used Eq. ([3]) with N = 10 6 for each set of param- 
eter values z, m and e, averaging over 10 3 random initial 
conditions xq uniformly distributed on the interval (0, 1), 
after discarding 4 x 10 3 iterates as transients. Figure 2 
shows the Lyapunov exponent for the map Eq. |I| as a 
function of the parameter e, for fixed z and several values 
of m, and for a fixed value of m and different values of 
the nonlinearity exponent z. 

The onset of type-Ill intermittency is associated to a 
direct transition from a fixed point to a chaotic one-band 
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FIG. 2: Lyapunov exponent A as a function of e for different 
values of m and with fixed z = 2; from top to bottom m = 
1, 1.5, 2, 3, 4, 6 and 8. Insert: A as a function of e, with fixed 
m — 2. for different values of z; from top to bottom z = 
2,3,4,5,6 and 8. 



attractor. The map Eq. ([T]) exhibits robust chaos for 
e > 0. A chaotic attractor is said to be robust if, for 
its parameter values, there exist a neighbourhood in the 
parameter space with absence of periodic windows (i.e. 
with A > 0) and the chaotic attractor is unique [13] • Ro- 
bustness is an important property in applications that 
require reliable operation under chaos, in the sense that 
the chaotic behaviour cannot be destroyed by arbitrarily 
small perturbations of the system parameters. It should 
be noted that robust chaos not associated to type-Ill 
intermittency has also been discovered in smooth, con- 
tinuous one-dimensional maps [l8|. 

The transition to chaos via type-Ill intermittency is 
manifested by a discontinuity of the derivative of the 
Lyapunov exponent at the bifurcation value e = 0. This 
discontinuity is due to the sudden loss of stability of the 
fixed point associated to the inverse period doubling bi- 
furcation that occurs at the onset of type-Ill intermit- 
tency. The Lyapunov exponent can be regarded as an or- 
der parameter that characterises the transition to chaos 
via type-Ill intermittency. This transition can be very 
abrupt as seen in [21 Thus, the behaviour of the Lya- 
punov exponent for e — > can be described by a scaling 
relation 



A - e 0{z > m \ 



(4) 



where /3 is a critical exponent expressing the order of 
the transition and that depends on the parameters z and 
m. Figure GJa) shows a log-log plot of the Lyapunov 
exponent vs. e for a fixed value of z and several values of 
the reinjection parameter m, while Fig. [3]Jb) shows the 
same plot with a fixed value of m and for different values 
of the nonlinearity z. Thus, relation Eq. (HJ) is satisfied 
independently for each parameter z and m. The critical 
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FIG. 3: (a) Log-log plot of the Lyapunov exponent A vs. e 
with fixed z — 2 for different values of m; from top to bottom 
m — 1,2,3,4,5,7 and 9. (b) Log-Log plot of A vs. e with 
fixed m = 1 for different values of z; from top to bottom 
z = 2,2.4, 3,4,6,8 and 10. 

exponent (3 as a function of z and as a function of m can 
be calculated from the slopes of each curve in Fig. 02a) 
and[3fb), respectively. 

Figure [4] shows the resulting graph of the critical expo- 
nent (3 as a function of both z and m. The values of (3 are 
distributed continuously on the interval (0, 1). We find 
the value (3 = 0.495±0.005 (typical error) for the parame- 
ter values z = 3 and m = 1, for which Eq. JT]) corresponds 
to the tangent map originally associated to type-Ill in- 
tcrmittency and having uniform reinjection probability, 
as assumed by Pomeau and Manneville [l|, |2|. However, 
for other values of z and m there is a manifested devia- 
tion from the Pomeau-Manneville's prediction (3 = 1/2. 



The calculated range for the critical exponent includes 
the particular value of j3 « 0.79 obtained for a perturbed 
logistic map showing type-Ill intermittency [13]. 




FIG. 4: Critical exponent /3 as a function of 2 and m. 



3. Conclusions 

We have characterised the scaling behaviour of the 
Lyapunov exponent at the transition to chaos in a generic 
map exhibiting type-Ill intermittency by means of a crit- 
ical exponent (3. We have shown that (3 depends on both, 
the order of the tangency of the map, described by the 
parameter z, and the probability of reinjection into the 
laminar region, expressed by the parameter m. We have 
calculated numerically the critical exponent on the space 
of parameters (z, m) and have obtained that (3 varies 
smoothly on the interval (0,1), in contrast with earlier 
predictions for the scaling behaviour of the Lyapunov ex- 
ponent in type-Ill intermittency. It is to be expected that 
the scaling behaviour shown here for the Lyapunov ex- 
ponent could also be detected in experimental situations 
with type-Ill intermittent chaos, where deviations from 
the scaling properties of the average laminar length have 
already been found. 
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